Opinion formation in the kinetic exchange models: Spontaneous symmetry breaking 

transition 
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We propose a minimal multi-agent model for the collective dynamics of opinion formation in the 
society, by modifying kinetic exchange dynamics studied in the context of income, money or wealth 
distributions in a society. This model has an intriguing spontaneous symmetry breaking transition 
to polarized opinion state starting from non-polarized opinion state. In order to analyze the model, 
we introduce an iterative map version of the model, which has very similar statistical characteristics. 
An approximate theoretical analysis of the numerical results are also given, based on the iterative 
map version. 
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I. INTRODUCTION 

Recently physicists have been studying social phenom- 
ena and dynamics leading to thegrowth of the interdis- 
ciplinary field of "Sociophysics" [l[ . One of the problems 
is of "opinion formation" , which is a collective dynam- 
ical phenomenon, and as such is closely related to the 
problems of competing cultures or languages 05 Si- It 
deals with a "measurable" response of the society to e.g., 
political issues, acceptances of innovations, etc. Numer- 
ous models of competing options have been introduced to 
study this phenomenon, e.g., the "voter" model (which 
has a binary opinion variable with the opinion alignment 
proceeding by a random choice of neighbors) Q, or the 
Sznajd-Weron discrete opinion formation model (where 
more than just a pair of spins is associated with the deci- 
sion making procedure) [5| . There have been other stud- 
ies of systems with more than just two possible opinions 
Q , or where the opinion of individuals is represented by 
a "continuous" variable [S-Q using real numbers. Also, 
since opinion formation in a human society is mediated 



by social interactions between individuals, such social dy- 
namics has been considered to take place on a network of 
relationships (see ^2] for recent review on such models). 
A two body exchange dynamics has already been de- 
veloped in the context of modelling income, money or 
wealth distributions in a society p1]hT4|. The general 
aim was to study a many-agent statistical model of closed 
econoniy (analogous to the kinetic theory model of ideal 
gases) ITsI, where TV agents exchange a quantity x, that 
may be defined as wealth. The states of agents are char- 
acterized by the wealth {xi}, i = 1,2,. . . ,N, such that 
Xi > 0, Vi and the total wealth W = J2i ^i is con- 
served. The question of interest is: "What is the equilib- 
rium distribution of wealth /(x), such that f{x)dx is the 
probability that in the steady state of the system, a ran- 
domly chosen agent will be found to have wealth between 
X and X + dx?" The evolution of the system is carried 
out according to a prescription, which defines the trad- 
ing rule between agents. The agents interact with each 
other through a pair-wise interaction characterized by a 
"saving" parameter A, with < A < 1. The dynamics of 
the model (CC) is as follows 15]: 
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x'i = Xxi + e(l - X){xi + Xj) , 

c'j = Xxj + (1 - e)(l - X}{xi + Xj) , 



(1) 



where e (0 < e < 1) is a stochastic variable, changing with 
time. It can be noticed that in this way, the quantity x 



is conserved during the single transactions: x\ + a;'- = 
Xi + x^ , where x[^ and x'- are the agent wealths after the 
transaction has taken place. In general, the functional 
form for steady state distribution /(x) is seen to be close 
to the F-distribution [l^ [13] • As a further generalization, 
the agents could be assigned different saving propensities 
and the steady state distribution j{x) show Pareto-like 
power- law behavior asymptotically [ISl ll9| . 

Earlier, Toscani [20| had introduced and discussed ki- 
netic models of (continuous) opinion formation involving 
both exchange of opinion between individual agents and 
diffusion of information. Based on this model, During et 
al 21] proposed another mathematical model for opinion 
formation in a society that is built of two groups, one 
group of ordinary people and one group of strong opin- 
ion leaders. Starting from microscopic interactions among 
individuals, they arrived at a macroscopic description of 
the opinion formation process. They discussed the steady 
states of the system, and extended it to incorporate emer- 
gence and decline of opinion leaders. 

Here, we report the studies of a minimal model for the 
collective dynamics of opinion formation in the society, 
based on kinetic exchanges. 



II. MODEL FOR OPINION FORMATION AND 
RESULTS 

A. Homogeneous multi-agent model 

Following the CC model described in the earlier sec- 
tion, we present a minimal model |22| for the collective 
dynamics of opinion Oi (t) of the i-th person in the soci- 
ety, consisting of N {N — > 00) persons. We assume that 
any particular person can discuss (interact) only with one 
other person each time (time increases discretely by unity 
after each such discussion). A two-person "discussion" 
is viewed here as a simple two-body scattering process 
in physics. Persons in the society may bump onto each 
other randomly and exchange opinions through such ran- 
dom two-person discussions. In general, a person i could 
have any opinion Oi between two extreme polarities de- 
noted by -|-1 and -1. In any discussion at time i -I- 1, a 
person retains a fraction of his/her older opinion Oi (t) , 
determined by his/her "conviction", parameterized by 
Xi. This parameter value is characteristic of a person 
and does not change with time t. Additionally, the per- 
son i is "influenced" stochastically by the other person 
j during the discussion having the "influence" parame- 
ter equal to his/her conviction parameter Xj . We further 
assume for simplicity that all agents are homogeneous - 
have the same conviction parameter A. Mathematically 
the dynamics may be represented by 




0,it + l) = X{0,{t) + etOj{t)), 
0,(t + l) = A(0,(i)+40,(i)), 



(2) 



where the opinion —1 < Oi{t) < 1 for all agents i and 
time t, the conviction parameter < A < 1 is quenched 



FIG. 1. Numerical results for the variation of the average 
opinion 0(t) for large t (steady state value of O) against A, 
following dynamics of Eq. ^. (Inset) Numerical results for 
the variation of the variance (AO)^ = (O ~ O)"^ against A, 
following dynamics of Eq. ((2]) . 



(does not change with time) , and the stochastic parame- 
ters et and e^ are annealed variables (change with time)- 
uncorrelated random numbers uniformly distributed be- 
tween zero and unity. Note that the equations are linear, 
but non-linearity is introduced in this model by imposing 
that — 1 < Oi(i) < 1 for all agents i and times t. 

The question we are interested is that if such social 
dynamics continually take place, can any consensus be 
reached or polarity evolve after a long time? Mathemat- 
ically, we are interested in the steady state distribution 
of O and other statistical properties. It is noteworthy 
that unlike in the market models, here we have no con- 
servation of opinion. Rather, the steady state of value of 
0{t) = {l/N)\YjiOi{t)\ represents the order of the aver- 
age opinion in the society after a long time t. We study 
the relaxation dynamics in the society: the relaxation 
and fluctuation of O, the steady state value of 0{t) for 
t > T, the relaxation time. 

Remarkably, we find there is appearance of "polar- 
ity" or consensus, starting from initial random disorder 
(where O^'s are uniformly distributed with positive and 
negative values). In the language of physics, there is a 
"spontaneous symmetry breaking" transition in the sys- 
tem: starting from (5(0) = the system evolves either 
to the "para" state with O = 0{t > r) = (where 
all individual agents have the opinion 0) for A < 2/3, 
or (continuously) to the "symmetry broken" state = 
0{t > t) ^ (where all individuals have either posi- 
tive or negative opinions) for A > 2/3 (see Fig. [Ij for 
times t > T. We note, however, that the fluctuation in 
O does not diverge, and shows a cusp near Ac (see inset 
of Fig. [1]). We also study the relaxation behavior of 0{t) 
and the critical divergence of the relaxation time r near 
A = Ac = 2/3 (see subsection HTCl Fig.H]). 
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FIG. 2. Numerical results for the variation of the average 
opinion 0{t) for large t (steady state value of O) against A, 
following dynamics of Eq. Q. (Inset) Numerical results for 
the variation of the variance (AO)^ = (O — O)^ against A, 
following dynamics of Eq. ^ . 



B. Random multiplier map 

The basic nature of transition produced by Eq. ^ , can 
perhaps be reproduced by the following simple iterative 
map 



Oit + I) = X{1 + et)0{t) 



(3) 



with the restriction that 0{t) < 1, which is ensured by 
assuming that if 0{t) > 1, 0{t) is set equal to 1. As 
usual, £( is a stochastic variable ranging between and 
1 (assumed to be uniformly distributed in our case). In 
a mean-field approach, the above equation reduces effec- 
tively to a multiplier map like 0{t -I- 1) = A(l -|- {e))0{t), 
where (e) = 1/2. Clearly for A < 2/3, 0{t) converges to 
zero. The initial value 0(0), is assigned either a positive 
or negative value. If it starts from a positive (negative) 
value, 0{t) remains positive (negative). We note that 
there are subtle differences in the dynamics of Eq. ([2|) and 
Eq. ([3]). Apart from the absence of "spontaneous sym- 
metry breaking" of the multi-agent model (from ±0^(0) 
values to all positive or all negative transition beyond 
Ac), the nature of the phase transition (singularity) in 
the iterative map is also slightly different. The critical 
value Ac = exp{ — (21n2 — I)} « 0.6796 has an analyti- 
cal derivation [23|, but for most numerical studies done 
here, we take Ac = 0.68. The time variation of the av- 
erage opinion 0{t) = {l/N)'Ei\Oi{t)\, where i refers to 
different initial realizations and N refers to the total of 
all such realizations, and its fluctuations are studied nu- 
merically (see Fig.[2|). We study the relaxation behavior 
of 0{t) and the critical divergence of the relaxation time 
T near A = Ac = 0.68 (see subsection III CI Fig. g]). We 
note again that the fluctuation in O does not diverge, 
and shows a cusp near Ac — 0.68 (see inset of Fig. [2]). We 
also note that the steady state fluctuation AO near Ac, 
is generally much higher in magnitude for the map case. 



C. Results and analyses 

For both the multi-agent model and the iterative map, 
we study the variation (with A) of the fraction p of the 
agents having O^ = ±1 at any time t in the steady state 
{t > r) . This parameter p gives the average "condensa- 
tion" fraction (of people in the society having extreme 
opinions \Oi\ = 1) in the steady state. The growth of 
p, as shown in Fig. [31 is seen to be similar to that of 
O. The inset shows that the growth behavior for p above 
(respective) Ac, for both the multi-agent model and map, 
are identical. 

We studied the relaxation behavior of O and p, for both 
the multi-agent model and map. In each case, the relax- 
ation time is estimated numerically from the time value 
at which O oi p first touches the steady state value within 
a pre-assigned error limit. We find diverging growth of 
relaxation time r (for both O and p) near A = Ac (see 
Fig. S]). The values of exponent z for the divergence in 
T ~ |A — Ac|~^ have been estimated numerically for both 
the multi-agent model and the map (for both A > Ac 
and A < Ac, wherever accurate data were obtained). For 
the multi agent model, the fitting values for exponent z 
corresponding to O and p, respectively, are z w 1.0 ± 0.1 
and z w 0.7 ± 0.1. For the map case, the fitting values 
for exponent z corresponding to both O and p, turn out 
to be the same: z « 1.5 ± 0.1. 

For the iterative map Eq. (jH]), we study carefully the 
time evolution of the condensation fraction p of |0| = 
1 in different realizations at different values of A. The 
variation of the steady state value p against A is shown 
in Fig. [21 It may be noted that while the steady state 
value of O starts to grow from A « 2/3 (see Fig. [2]) , the 
steady state value of p starts growing at A « 0.68 (see 
Fig. [3]). Numerical results for the growth of the relaxation 
time r for both O and p, against A are shown in Fig. [J 
Both diverge at A w 0.68. This clearly indicates that p, 
rather than O, is the order parameter for the transition. 

An approximate analysis of the above transition for A 
closer to unity can be done for the iterative map Eq. ([3]) 
as follows. In Fig.[5l we give the numerical results for the 
steady state distribution opinion, P(|0|) for three differ- 
ent values of A; we observe roughly a bi-modal nature 
of the distribution as A — >■ 1: one mode is the uniform 
distribution within the range \Omin\ < \0\ < 1 (and 
\Omm\ ~ A) and another a delta function at \0\ — 1. 
We therefore approximate the steady state distribution 
of opinion by assuming that opinion 0{t) is distributed 
uniformly starting from a minimum Omin upto unity with 
(integrated) probability (1 ~p), and a J- function at ex- 
actly unity with probability p. Then 

O = {1 - p)Oa. + p.l, (4) 

where Oav = {Omin + l)/2. We have assumed that the 
value 0{t) stays in those two regions (from A to 1 and 
exactly at 1) with probability (1 — p) and p. Hence, the 
corresponding equations are 



0(i -t- 1) = A(l + e)O(i) 



with probability 1 ~ p, 
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FIG. 3. Numerical results for the variation of the average 
condensate fraction p(t) for large t (steady state value of p) 
against A, following dynamics of Eq. ([2]) in black diamonds, 
and dynamics of Eq. Q in red circles. (Inset) Numerical 
results for the growth of p, following dynamics of Eq. ((2} in 
black diamonds, and dynamics of Eq. ^ in red circles, close 
to Ac. 



and 



0(i + l) = l 



with probability p. 



^Oav 



Note that the first equation is realized only if A(l 

e)O(i) < 1 or e < emax ~ j^ 1- This cut-ofF implies 

that p = J^""" 
substituting Oa 
that 



de 



1, since e '^uni[0,l]. By 
and p in Eq. (J4l), we derive the result 



AOa„ 



= 



5A + 2A2 - A3 - 2 
2A(H-A) 



(5) 



which is compared with the numerical simulations for 
A ^^ 1 in Fig. |6l It is evident that the approximation 
holds well, only for A — s> 1. 



III. SUMMARY AND DISCUSSION 

In summary, we proposed a minimal model for the 
collective dynamics of opinion formation in the society, 
by modifying kinetic exchange dynamics studied in the 
context of markets. The multi- agent model (dynamics 
given by Eq. ([2])) and its map version (dynamics given 
by Eq. (|31)) have kinetic exchange like linear contribu- 
tions from random two-person discussions or scattering 
processes, though the saturation of \Oi\ < 1 induces non- 
linearity in the dynamics. This model has an intriguing 
spontaneous symmetry breaking transition to polarized 
opinion state starting from non-polarized opinion state. 
Specifically, in the multi-agent model, we see that for 
A > Ac = 2/3, starting from random positive and neg- 
ative Oi values (or for that matter any arbitrary state), 
at t = 0, the system eventually evolves to a state at 
t > T where all Oi are either positive or negative, with 
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\0\ determined by the A value! This is similar to the 
growth of spontaneous magnetization in Ising magnets 
(where starting from arbitrary up and down spin states, 
a preferred direction is chosen by fluctuation) , with mag- 
netization determined by the temperature below its tran- 
sition value. The appearance of spontaneous symmetry 
breaking in this simple kinetic opinion exchange model 
is truly remarkable. It appears to be one of the sim- 
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FIG. 5. Numerical results for for the steady state distribution 
opinion, P{\0\), for three values A = 0.8, 0.85, 0.9 showing bi- 
modal distributions in each case. (Inset). The same steady 
state distribution PflOl). for three values A = 0.8.0.85.0.9. 
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plest collective dynamical model of many-body dynamics 
showing non-trivial phase transition behaviour. Indeed, 
it may be noted that for A < Ac, at t > t, all Oi's 
become identically zero (without any fluctuation), while 
for A > Ac, Oi's have fluctuations but the average has a 
steady state value depending on the value of A. As such 
the nature of the phase transition in this model is quite 
different and does not fit to the commonly studied three 
absorbing state models studied (see e.g., [2J|)- In order to 
understand the nature of the transition, we also studied 
a simple iterative map and derived approximate result 
for the order parameter variation under certain limits, 
which compares quite well with the numerical simula- 
tions. Specifically, we find that the fraction p of people 
with extreme opinion \0i\ — 1, and its fiuctuations de- 
termine the nature of the phase transition in our model 
and locate the critical point accurately (from numerical 
studies). With the two mode distribution (uniform and 
delta) of O, valid close to A ^^ 1 (see Fig. [5]), we could 
develop an approximate analysis of the variation of the 
steady state mean opinion \0\ against A as in Eq. ^. 
In any case, further investigations are necessary for un- 
derstanding this phase transition. Additional studies for 
the heterogeneous conviction factors A^'s, in influence of 
"field terms" that represent the external influence of me- 
dia, etc. will be reported elsewhere [25]. Also, the study 
of this phase transition behavior for an extended model 
with separate "conviction" and "influence" parameters in 
Eq. ^, has recently been reported [2a |. 
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FIG. 6. Fit of the approximate theoretical calculation Eq. ([5]) 
(in red squares) with the numerical simulations for A — >■ 1, 
following dynamics of Eq. Q (in black circles) . 
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